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Abstract 

We derive the evolution of the hnear bias factor, h{z), in cosmological models driven 
by an exotic fluid with an equation of state: px = wpx, where — 1 < w < (quintessence). 
Our aim is to put constrains on different cosmological and biasing models by combining the 
recent observational clustering results of optical {2dF) galaxies (Hawkings et al.) with those 
predicted by the models. We find that our bias model when fitted to the 2dF clustering 
results predicts different bias evolution for different values of w. The models that provide 
the weak biasing {ho ~ 1.1) of optical galaxies found in many recent observational studies 
are fiat, Vt^ = 0.3 with w < —0.9. These models however, predict a weak redshift evolution 
of b{z), not corroborated by N-body simulations. 
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1. Introduction 

The statistical properties of the distribution of matter on large scales, based on differ- 
ent extragalactic objects, can provide important constrains on models of cosmic structure 
formation. However, a serious problem that hampers the straight forward use of such an 
approach is our ignorance of how luminous matter trace the underlying mass distribution. 
Many authors have claimed that the large scale clustering of different mass tracers (galaxies 
or clusters) is biased with respect to the matter distribution (cf. Kaiser 1984), which is an 
essential ingredient for cold dark matter (CDM) models to reproduce the observed galaxy 
distribution (cf. Davis et al. 1985). Usually, biasing is assumed to be statistical in nature by 
which galaxies and clusters are identified as high peaks of an underlying, initially Gaussian, 
random density field. 

Furthermore, the bias redshift evolution is very important in order to relate observations 
with models of structure formation and it has been shown that the bias factor, b{z), is a 
monotonically increasing function of redshift. Indeed, Mo & White (1996) and Matarrese et 
al. (1997) have developed a model for the evolution of the correlation bias, defined as the ratio 
of the halo to the mass correlation function, the so called galaxy merging bias model, with 
b{z) cx (1 + zy-^. Nusser & Davis (1994), Fry (1996), Tegmark & Peebles (1998) produced a 
bias evolution model assuming only that the test particles fluctuation fleld is proportional to 
that of the underlying mass and in this case, we have: b{z) — 1 + (6o — 1)(1 + where bo 
is the bias factor at the present time. Coles, Melott & Munshi (1999) have developed a bias 
model within the hierarchical clustering paradigm which gives interesting scaling relations 
for the galaxy bias. Recently, another approach was proposed by Basilakos & Plionis (2001), 
in which the bias evolution was described via the solution of a second order differential 
equation, derived using linear perturbation theory and the basic Cosmological equations. 

From the observational point of view, recent advances in cosmology have strongly indi- 
cated that we are living in a flat, accelerating universe with low matter (baryonic and dark 
matter) density (cf. Riess et al. 1998; Perlmutter et al. 1999; dc Bernadis et al. 2000; 
Efstathiou et al. 2002; Persival et al. 2002, Spiergel et al 2003). The available high quality 
cosmological data (Type la supernovae, CMB, etc.) are well fltted by an emerging "standard 
model", which contains cold dark matter (CDM), to explain clustering, and an extra com- 
ponent with negative pressure, usually named "quintessence", which is in agreement with 
the inflationary flatness prediction (Jltot = 1) as well. 

In this paper we derive the redshift evolution of the linear scale-independent bias pa- 
rameter for the "quintessence" models, and additionally we attempt to put constraints on 
the different cosmological models using the recent clustering results (Hawkins et al. 2002) 
of optical galaxies from the 2dF survey. We compare the observed clustering of 2dF galaxies 
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with that expected in 10 spatially flat cosmological models with dark energy and = 0.3. 
It is worth noting that the galaxy bias has been found to be sensitive to the equation of state 
(Munshi, Porciani & Wang 2003). 



In this section we briefly describe our linear bias evolution model (for more details 
see Basilakos & Plionis 2001) which is based on hnear perturbation theory in the matter 
dominated epoch (cf. Peebles 1993). We assume that the mass tracer population is conserved 
in time, ie., that the effects of non-linear gravity and hydrodynamics (merging, feedback 
mechanisms etc) can be ignored (cf. Fry 1996; Tegmark & Peebles 1998; Catelan et al. 
1998), and using the linear perturbation theory we can obtain a second order differential 
equation, which describes the evolution of the linear bias factor, b, between the background 
matter and a mass-tracer fluctuation fleld: 



where D{t) is the linear growing mode, useful expressions of which can be found for the 
ACDM model in Peebles et al. (1993) and for the QCDM models in Wang & Steinhardt 
(1998). The solutions of this second order differential equation provides our bias evolution 
model. The solution for the different cosmological model enter through the different be- 
haviour of D{t) and H{t). To see the relevance and difference of our bias evolution model to 
the Galaxy Conserving one, which is based on similar assumptions, we will derive somewhat 
differently eq.(l), starting from the continuity equation, which is the starting point of the 
Galaxy Gonserving model. If the galaxies and the underlying matter share the same velocity 
field then: 6 + Vm ~ and 6g + Vu ~ 0, from which we get d — Sg = 0. Now since we assume 
linear biasing, we have 6g = b6 and using y = b—1, v/e get that d{y6)/dt = 0. Differentiating 
twice we then have: y6 + 2yS + y5 — 0. Solving for y5, using the fact that y5 — —y5 and 
utilizing the differential time-evolution equation of 6 (cf. Peebles 1993) we finally obtain: 



which is the corresponding equation 1. 

Here we extend our previous bias evolution solutions of Basilakos & Plionis (2001) to 
flat Friedmann-Robertson Walker Cold Dark Matter (CDM) type models driven by non- 
relativistic matter and having an exotic fluid (quintessence) with equation of state: 
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Linear Bias Evolution Model 



bD{t) + 2 D{t) + H{t)D{t) b + AnGp^D{t)b = inGpraD{t) 



(1) 



yS + 2{S + HS)y + AirGpSy = 



(2) 



Px = WPx (-1 <W <0) 



(3) 
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FoUowing the notation of Basilakos & Plionis (2001), we present the basic steps of our 
procedure. In order to transform (1) from time to redshift we utihze the following expressions: 

I = -imkri) • ^(^' - ["-<^ + + "-<^ + ^'^ 

while Qm = ^TiGpo/'iH^ (density parameter), Q^^ = STrGp^^/SiJ^ (dark energy parameter) 
at the present time, which satisfy + Q^o = 1; /9 = 3(1 + w) with < /9 < 3. Taking 
into account the latter transformations, the basic differential equation for the evolution of 
the linear bias parameter takes the following form: 

d% .d6 



with relevant factors, 



^ ' 1 + z E{z) dz D{z) dz ^ ^ 

and 

Q(z) = ^^-(^ + ") (7) 

Therefore, the general bias solution for all of the flat cosmological models is: 

b{z) = AE{z) + CE{z)I{z) + 1 , (8) 

where 

The integral of equation (9) is elliptic and therefore its solution, in the redshift range [z, +oo), 
can be expressed as a hyper-geometric function. Our final general solution is: 

1 3 7-2/5 



2C 13-3 

biz) - 1 = AE{z) + ^E{z){l + z)^F 



(10) 



_6-2/?'2' 6-2/?' 0^(l + z)3-/3_ 

Note that the first term, which is the dominant one, has an approximate redshift dependence 
~ (1 + z)^/^ while the second has ~ {1 + z). Also note that for w — > —1 the above general 
bias solution tends to the ACDM it should (eq.39 of Basilakos & Plionis 2001). It is 

quite instructive to see the solution for the case of tt; = — 1/3, in which case the functional 
form of the bias evolution is independent of the hyper-geometric function^. We have that: 

b{z)^Ail + z){l + n^zy/^ + -^{l + z) + l . (11) 

We can easily verify that for Q^a — 1 this reduces to the Einstein-de Sitter case (eq. 29 of 
Basilakos & Phonis 2001), as it should. 



^Where we have used 0^; = 1 — fl^n and F[a, b, b, t] = {1 — t) 
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3. Estimating the Bias from Galctxy and Mass correlations 



Since our approach gives a family of bias curves, due to the fact that it has two unknown 
parameters, (the integration constants A, C) and in order to obtain partial solutions for b{z) 
we need to estimate the values of these constants. In Basilakos & Plionis (2001) we have 
compared our bias evolution model with the halo merging model (cf. Mo & White 1996; 
Mataresse et al. 1997) as well as with different N-body results and found a very good 
consistency, once we fitted the integration constants A,C by evaluating our model to two 
different epochs. We further compare in figure 1 our solution for the ACDM case (see its 
parameters further below), evaluated at z = and z = 3 using the HDF results (Arnouts et 
al. 2002; Malioccietti 1999, with the Mataresse et al (1997) model. It is quite evident that 
our model fits better the 2;-dependance of the observational HDF galaxy bias. 

Although, this comparison gives consistent results of the functional form of our solution 
with available data and theoretical models (see Basilakos & Plionis 2001), it does not test 
directly whether our model, once calibrated observationally, provides a consistent model at 
high z's. Therefore, in order to test it, we evaluate the constants A,C by fitting our bias 
model to the the recent 2dF galaxy clustering results of Hawkins et al. (2002) who used 
~200000 galaxies observed in the 2dF survey to derive their spatial correlation function. 

We use the standard theoretical approach to estimate the two point spatial correlation 
function, using our model for the bias evolution in the different spatially flat cosmological 
models. We quantify the evolution of clustering with epoch, writing the spatial galaxy cor- 
relation function as ^modeK'^; 2;) = Cmass('")-R(^); with R{z) = D^{z)b'^{z), where the function 
R{z) characterizes the clustering evolution with epoch (see Basilakos 2001 and references 
therein). While the ^mass('^) is the Fourier transform of the spatial power spectrum P{k): 



where k is the comoving wavenumber. Note that we also use the non-linear corrections 
introduced by Peacock & Dodds (1994). As for the power spectrum of our CDM models, 
we take P{k) ^ k^T'^{k) with scale-invariant {n — 1) primeval inflationary fluctuations and 
T{k) the CDM transfer function. In particular, we use the transfer function parameterization 
as in Bardeen et al. (1986), with the corrections given approximately by Sugiyama's formula 
(Sugiyama 1995). 

In the present analysis we consider flat models with cosmological parameters that flt 
the majority of observations, ie., + — 1; = 0.3, Ho — 100/ikm s~^ Mpc~^ with 
/i ~ 0.6 — 0.7 (Freedman et al. 2001; Peebles and Ratra 2002 and references therein) and 
baryonic density parameter ^5/1^ ~ 0.02 (cf. Olive, Steigman & Walker 2000; Kirkman et 
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al 2003). We investigate such spatially flat cosmological models with negative pressure for 
different values of w (see eq. 3). To this end, all the cosmological models are normalized to 
have fluctuation amphtude in 8 h~^Mpc scale of erg = 0.50(±0.1)O~''' (Wang & Steinhardt 
1998) with 7 = 0.21 - 0.22w + 0.33fl^. 

In order to quantify the integration constants {A,C) we perform a standard mini- 
mization procedure between the measured correlation function for the 2dF galaxies (Hawkins 
et al. 2002) with those expected in our spatially flat cosmological models, 



where a is the observed correlation function uncertainty. 

In figure 2 we present the fit to the data of two models with values of w at the opposite 
end of the range explored. Hawkins et al (2002) found that the 2dF clustering behaviour 
on small scales is represented by a power law, with correlation length Tq ~ 5.05/i~^Mpc 
and slope 7 ~ 1.67. Due to the interplay between the three unknown {A,C,bo) all of the 
models (for different values of the constants) fit at a high significance level (P^2 ^0.7) the 
2dF correlation function. In Table 1 we list the results of the fits for all models, ie., the 
integration constants and the value of the optical bias, 60, at the present time, derived from 
our model (eq. 10). It is interesting that for the standard ACDM model {w — —1, Qrn — 0.3 
and h — 0.7) the ~ C{l + z) factor of the solution is ~ 0, which means that the bias evolution 
is a very weak function of z, approximated by: b{z) — 1 ^ ^(1 + z)^/^. 

It is also very interesting that for the ACDM and the w — —0.9 QCDM models (with 
h — 0.7), the corresponding bias at the present time is bo ^ 1.05 and 60 — 1-11, respectively, 
in very good agreement with the values bo — 1.04 and 1.11 derived by Verde et al. (2002) 
and Lahav et al. (2002), respectively. When we use h = 0.6. only the ACDM model provides 
a value for the present bias which is only roughly consistent with other observational results 
(bo — 1.17). The present time bias, bo, as a function of w, is well fitted by a power law 
having the form: bo ~ [1.93 - 1.20/i] (±0.019) x (_^i;)-o-4ii(±o.oi8)_ 

Having derived the values of the constants A,C and the zero-point bias, 60, from the 
fit to the 2dF galaxy clustering pattern, we can now investigate the bias evolution in the 
different cosmological models. In figure 3 we present the function b{z) for different values 
of w. The biasing is a monotonically increasing function of redshift with its evolution being 
significantly stronger for lower values of w. Galaxy clustering in quintessence models evolves 
more rapidly than in the ACDM model. It is also clear that for the preferred cosmological 
models (with w < —0.9 and h — 0.7) the optical galaxies are only weakly biased, even at 
high redshifts. 
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In figure 3 we also compare our evolution model with the galaxy conserving one (Pry 
1996; Tegmark & Peebles 1998) by normalizing the latter to the value of ho provided by our 
fit. It is interesting that the two models, despite their different form, are consistent with 
each other. However, we note that this is true for the h = 0.7 case. Had we presented the 
h = 0.6 case it would have be evident that the consistency is significantly worse. In the same 
figure we compare our model with the results of a ACDM N-body simulation of Somerville et 
al (2001) (thick blue line). The simulation galaxy bias evolves significantly more than what 
our model predicts, which should be attributed to our assumption that the galaxy number 
density is conserved in time. It is evident that merging processes, not taken into account by 
our model, are very important in the evolution of clustering. 



4. Conclusions 

In this paper we derived analytically the evolution of linear bias b{z) for cosmological 
models driven by an exotic fluid with an equation of state: — wp^, where —l<w<0. 
Comparing the optical 2dF galaxy correlation function (Hawkings et al. 2002) with the 
predictions of various (quintessence) models, we find that the fiat cosmological models with 
w < —0.9 provide present bias values in the range 1.05 < bo < 1.11, which is consistent with 
observational correlation results. However, we find that for the above cosmological models 
our bias evolution is too weak to fit N-body results of galaxy clustering, a fact which should 
be attributed to merging processes, not taken into account in our model. 
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Fig. 1. — Comparison of our bias model for the ACDM case (continuous line) with the 
corresponding Mataresse et al (1997) halo merging model for galaxies with ~ 10^^ Mq 
(dashed line). The points represent the HDF galaxy bias data from the HDF- north (open 
symboles) and HDF-south (filled symboles). Note that our model has been normalized to 
the data at ^ = and z = 3. 
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Fig. 2. — Comparison of the 2dF correlation function (Hawkins et al. 2002) with that 
predicted for the ACDM {w = —1) and an extreme quintessence model {w = —0.1), using 
h — 0.7. Each fit corresponds to different values of the integration constant A, C and of the 
zero-point bias, bo (see table 1). 
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Fig. 3. — Evolution of the linear bias factor for our model with h = 0.7. Different line types 
represent different values of w in our model, while the filled points is the Fry (1996) evolution 
model normalized to b{0). The thick continuous line (blue) represent the results of galaxies 
with Mb — 5\ogh < —19.5 in the Somerville et al (2001) ACDM simulation while the star 
like symboles represent the HDF results for the ACDM model of Maliochetti & Lahav (1998). 
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Table 1: Our bias model normalizations and present linear bias {bo) between optical galaxies 
and underlying matter distribution. 





Ho 


= 70 km s-^Mpc" 


-1 


Ho 


— 60 km s ^Mpc" 


-1 
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A 
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bo 


A 


C 


bo 


-1.0 


0.016 


0.004 


1.05 


0.004 


0.020 


1.17 


-0.9 


0.028 


0.010 


1.11 


0.022 


0.028 


1.23 


-0.8 


0.038 


0.018 


1.18 


0.014 


0.038 


1.31 


-0.7 


0.024 


0.032 


1.27 


0.032 


0.050 


1.41 


-0.6 


0.084 


0.040 


1.37 


0.054 


0.064 


1.53 


-0.5 


0.074 


0.064 


1.50 


0.058 


0.092 


1.67 


-0.4 


0.190 


0.078 


1.67 


0.070 


0.130 


1.85 


-0.3 


0.190 


0.130 


1.86 


0.090 


0.190 


2.07 


-0.2 


0.490 


0.150 


2.10 


0.320 


0.250 


2.34 


-0.1 


0.480 


0.320 


2.39 


0.410 


0.440 


2.67 



